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ABSTRACT 

We present details of the construction and characterization of the coaddition of the Sloan Digital 
Sky Survey Stripe 82 ugriz imaging data. This survey consists of 275 deg 2 of repeated scanning by the 
SDSS camera of 2.5° of 8 over —50° < a < 60° centered on the Celestial Equator. Each piece of sky 
has ~ 20 runs contributing and thus reaches ~ 2 magnitudes fainter than the SDSS single pass data, 
i.e. to r ~ 23.5 for galaxies. We discuss the image processing of the coaddition, the modeling of the 
PSF, the calibration, and the production of standard SDSS catalogs. The data have r-band median 
seeing of 1.1", and are calibrated to < 1%. Star color-color, number counts, and psf size vs modelled 
size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure 
in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars 
are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide 
area deep survey there are a variety of uses for the data, including galactic structure, photometric 
redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass 
calibrations, and cosmic shear measurements. 
Subject headings: atlases — catalogs — surveys 



1. INTRODUCTION 



The Sloan Digital Sky Survey (SDSS; [YorkeFaL 2000) 
saw first light in 1998 with the goal or obtaining CCD 
imaging in five broad bands ugriz over 10,000 deg 2 of 
high-latitude sky in the North Galactic Cap, plus spec- 
troscopy of one million galaxies and one hundred thou- 
sand quasars over this same region. In addition, the 
SDSS imaged a 275 deg 2 region on the Celestial Equa- 
tor in the Southern Galactic Cap. This region is called 
"Stripe 82" and was imaged multiple times during the 
Fall months when the North Galactic Cap was not ob- 
servable. The SDSS single pass data reach r ~ 22.4 and 
has median seeing of 1.4" in r, but by aligning and aver- 
aging ( "coadding" ) the Stripe 82 images we reached ~ 2 
magnitudes deeper and median seeing of ~ 1.1". The in- 
tent was to use this deep survey to understand the single 
pass data at its limits and to do science at fainter magni- 
tudes or correspondingly higher redshifts. Such analyses 
benefit from our image processing approach, as opposed 
to catalog-level methods, because objects below the de- 
tection limit of individual single pass images can be de- 
tected and measured. A brief description of the Stripe 82 
data and coadd was pr esented in the SPSS Sev enth Data 
Release (DR7) paper (Abazajian et al. (2009); see also, 
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Jiang et al. (2008)). Here we give a full report, detailing 
the features in he coaddition process. 

The SPSS uses a dedicated wide-field 2.5m telescope 
(Gunn et al. 20061 located at the Apache Point Obser- 
vatory (APO) near Sacramento Peak in So uthern New 
Mexi co. The telescope imaging instrument ( Gunn et al 



1998) is a wide-field camera with 24 2048 x 2048 0.396" 
pixel scale CCPs. SPSS images the sky in drift scan 
mode with the five filters in the order riuzg (Fukugita 
et al. 1996 ) . Imaging is performed with the telescope 



tracking great circles at the sidereal rate; the effective 
exposure time per filter is 54.1 seconds, and 18.75 deg 2 
are imaged per hour in each filter. The images are 
mostly taken under goo d seeing condition s on moon- 



less photometric nights (Hogg et al. 20011. For stel- 
lar sources the 50% completeness limits of the images 
are u,g,r,i,z = 22.5, 23.2, 22.6, 21.9, 20.8, respectively 



( Abazajian et al.|2003 ), although these values depend on 
seeing and sky brightness. The image p rocessing pipelin e 
determines the astrometric calibration ( Pier et aL|2 003), 
then detects objects and measures their brightnes ses, po- 
sition s and shapes ( Lupton et al.||200T Stoughton et al. 
2002 1 . The astrometry is good to 45 milli arcseconds 
(mas) rms p er coordinate at the bright end (Abazajian 
et al.| |2 009[). The photome try is calibrated to an AH 



systern fOke fc Gunn||1983l, and the zero-points of th e 
system are known to 1-2% ( [Abazajian et al.|20 03 2004 ) 
The photometric calibration is done in two ways, by ty 



ing to photometric standard stars (Smith et al. 2002] ) 
measured by a separate 0.5m telescope on si te (the FT 
telescope; Tucker et al. 2006 Ivezic et al. 2004 ) and by us- 
ing the overlap between adjacent imaging runs to tie the 
photometry of all the imaging obs ervations together, in 
a pro cess called "ubercalibration" (Padmanabhan et al. 
20081. Ubercalibration zero-points on each stripe have 
rms error of ~ 2% in u and ~ 1% in griz. 

SPSS data is obtained as runs, where a run is a single 



2 



Annis et al. 



continuous drift scan obtained on a single night. A sur- 
vey stripe is one camera width wide, about 2.5°. Two 
interleaving runs called strips are necessary to complete 
a stripe as the camera focal plane is sparsely populated. 
These strips are denoted either N or S, depending if the 
telescope boresight is pointed half a CCD width north or 
south of the stripe equator. A run contains 6 columns 
of data through the five ugriz filters, and a single fil- 
ter data set is called a scanline. Each scanline is a 13' 
wide continuous stream of data that we arbitrarily chop 
into overlapping 10' long frames. A frame is a single im- 
age in a single bandpass, and has a geometry of 1489 
rows and 2048 columns, at a pixel scale of 0.396" /pixel. 
A field is the set of ugriz frames of the same piece of 
sky, disregarding the fact that they were obtained over 
8 minutes of time. For Stripe 82 in particular, there is a 
unique mapping of RA and Dec into survey constructs. 
The row number of a field corresponds to RA, as does 
the field number. The column number of a field corre- 
sponds to Dec, as does the camera column number. The 
fields overlap along the RA direction by 124 rows per 
field due to a repackaging of the same pixel data during 
data acquisition. They overlap along the Dec direction 
by a small amount on either edge of the field due to 
re-observation of sky by slightly overlapping scanlines. 
The coadd runs are artificial, so we adopted run number 
100006 as the south strip and 200006 as the north strip 
arbitrarily. These were later renamed as 106 and 206 for 
convenience. These runs interleave, but for constant col- 
umn number, run 206 is at a higher Dec than run 106. 
Each run is 800 fields long and goes in increasing field 
number from the West to East, from low to high RA. For 
more informati on on the SPSS nomen clature and tech- 
nical terms, see Stoughton et al.| ( 2002 1. 

Stripe 82 is the SDSS stripe along the Celestial Equa- 
tor in the Southern Galactic Cap. It is 2.5° wide and 
covers -50° < RA < +60°, so its total area is 275 deg 2 . 
Stripe 82 can be observed from APO at low airmass from 
September through November, is accessible from almost 
all ground-based telescopes for subsequent spectroscopic 
and photometric observations and , except near its RA 
ends, has low Galactic extinction (Schlegel et al. 1998). 
Stripe 82 was imaged by the SDSS multiple times in the 
Fall months and through 2004, these data were taken 
only under optimal seeing, sky brightness, and photomet- 
ric conditions (i.e., the c onditions required for imaging in 
the main Legacy Survey; York et al. (2000)). There were 
84 such runs. In 2005-2007, 219 additional imaging runs 
were ta ken on Stripe 82 as p art of the SDSS supernova 
survey ( Frieman et al.|[2008 ), designed to discover Type 
la supernovae at 0.1 < z < 0.4. The supernova survey 
was carried out on most usable nights, with the excep- 
tion of the five brightest nights around each full moon. 
Therefore, these data were often taken under less opti- 
mal conditions: poor seeing, bright moonlight, and/or 
non-photometric skies. 

Both reduced images and catalogs from all 303 runs 
covering St ripe 82 were made ava ilable as part of the 
SDSS DR7 ( Abazajian etaLl 2009|, in a database called 
Stripe 82. The data can be accessed both the Data 
Archive Server (DAS) and the Catalog Archive Server 
(CAS). We carried out a coaddition of the repeat imaging 
scans, photometric or not, on Stripe 82 taken through 
Fall 2005. Data taken after that date were excluded as 



for the most part they had not been taken when we were 
processing the coadd. The coaddition includes a total of 
123 runs, covering any given piece of the 275 deg 2 area 
between 20 and 40 times. The S and N strip runs are 
designated 100006 and 200006, respectively, in the DAS, 
and 106 and 206 in the CAS database. 

We designed the coaddition program so that the out- 
put image format allowed us to run the SPSS standard 
measurement code, PHOTO (Lupton et al.|2001||Stoughton 



et al. 2002 Lupton et al. |2012p , on the coadd images 



I'his was important because: a) PHOTO has algorithms 
which had been extensively tested by the SPSS collabo- 
ration over the years; and b) the resulting data products 
are conveniently structured for joint analyses and com- 
parisons with the single pass data. Our method considers 
the repeat scans of Stripe 82 to be noisy, distorted real- 
izations of the true sky. The aim is to make our best 
estimate of the true sky as it would have been seen by a 
perfect SPSS camera on a larger telescope. Starting with 
the list of runs on Stripe 82 taken from the start of the 
survey to the Fall 2005 season, those fields of reasonable 
seeing (FWHM), transparency (T), and sky noise (<r s ) 
were selected for use in the coadd. The individual runs 
were remapped onto a uniform astrometric coordinate 
system. Interpolated pixels (due, e.g., to cosmic rays or 
bad columns) in each individual run were masked and 
the sky was subtracted from each frame. The images are 
coadded with weights that depend on FWHM, T and a s , 
providing optimal signal ratio to noise for point sources. 
PHOTO relics on an accurate point spread functio n (PSF) 
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suitably weighted sum of the PSFs made by PHOTO for 
each run. The coadded images were run through PHOTO 
yielding the catalog made available in the CAS Stripe82 
database. 

When using the coadd data for science it is impor- 
tant, just as with the main survey, to use the various 
processing flags associated with each detected object to 
reject spurious objects and to select objects wit h reliable 
photometry (as recommended, for example, by Richards 
et al. (2002 1). Since the coadd data was run through the 
SPSS pipelines, the standard flag set is available for all 
objects. However, some objects at magnitudes < 15.5 
that are saturated do not have the saturated flag set, so 
we recommend a magnitude cut to avoid them. 

In this paper we first consider the observations (?J2]), 
then describe the coadd image (^3| and catalog creation 
(Q. In section <|5 we present the results of our quality 
assurance tests and explore the features of this data set, 
highlighting the improvements in depth and seeing due 
to the coaddition process. We discuss the applications 
and science outcomes in ^6] and conclude in <|7] 

2. OBSERVATIONS 

Figure [l] shows the number of observations as a func- 
tion of RA for runs 106 (S strip) and 206 (N strip) sep- 
arately. The total number of images reaches ~ 100 for 
the S strip (blue, top curve) and ~ 80 for the N strip 
(black, top curve). About 30% of those runs are cali- 
brated (red and black, bottom curves) in the sense that 
the infrared sky camera indicated a minimum of clouds 
and a extinction solution was obtained for the 20" PT 
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telescope data taken at the same night of the run. The 
final number of images used in the coadd, shown as thick 
green (N strip) and red (S strip) lines, varies from 15 to 
34. The selec tion criteria to achieve this final sample is 
described in §2.1| Based on the number of observations 
used, we expect the coadd to be ~ 2 mag deeper than the 
single pass data and to show a difference of ~ 0.4 mag 
in depth between the shallowest and the deepest regions, 
assuming that the signal-to-noise ratio increases as y/~N. 



100 




RA (degrees) 

Fig. 1. — RA distribution of Stripe 82 observations for both 
runs 106 and 206, corresponding to S and N strips respectively. 
The total number of images reaches ~ 100 for the S strip (blue, 
top curve) and ~ 80 for the N strip (black, top curve). Nearly 
30% of those are calibrated (red and black, bottom curves). The 
number of images selected for the coadd in both N (thick green) 
and S (thick red) strips varies from 15 to 34. We therefore expect 
to see a difference of ~ 0.4 mag in depth between the shallowest 
and the deepest regions of the coadd, assuming that the signal-to- 
noise ratio increases as \/N. See text for details on the selection 
criteria. 



2.1. Field Selection Criteria 

The data selection criteria are listed in Table □ We 
chose all runs on Stripe 82 with 125 < run < 5924, i.e., all 
data obtained in or before December 1 2005, demanding 
that the runs were either on the N or S strip and reject- 
ing ~ 10 runs that were not offset (strip labeled "O"), 
plus one run that is a crossing scan at ~ 45° inclination. 
We then select the fields in r band, requiring seeing bet- 
ter than 2", sky brightness less than 19.5 mag/arcsec 2 
and less than 0.2 mag of extinction. The sky bright- 
ness cut corresponds to 2.5 times the median sky of 150 
DN, allowing at most 0.5 mag increase in sky noise. The 
majority of the data is uncalibrated so we also require 
that the field has enough stars for our relative calibra- 
tion method to work. We cut on the r-band parameters, 
but rejected all the corresponding data in ugiz. This 
choice maximizes homogeneity across filters, though not 
optimization for a given filter. 

The fraction of fields passing the seeing and sky noise 



cuts are 91% and 84%, respectively, and 77.5% pass both 
cuts jointly. However, 100% of the fields in the stan- 
dard SDSS runs and 95% supernova fields pass the trans- 
parency cut and the demand to have enough stars to 
calculate the photometric scaling. This indicates that 
the SDSS had a high threshold for classifying a night 
as photometric. Overall, 1,124,075 frames were included 
in various fields of the coadd. The thick red and green 
lines in Fig. [T] shows the RA distribution of these selected 
frames. 

Table [2] summarizes the 123 runs included in the coadd 
from both main survey and supernova runs. 69 of these 
were calibrated (phot=l), in the sense that the infrared 
sky camera indicated a minimum of clouds and a extinc- 
tion solution was obtained for the PT telescope data. 
The remaining 54 runs were uncalibrated (phot=0). 
These were on average twice as long as the calibrated 
ones, as they were taken by the supernova survey which 
was unconcerned with photometricity or sky brightness. 

Not all of the images in the available overlapping runs 
were used in a given coadded frame. In order to prevent 
sharp discontinuities in the PSF in the output image at 
input run edges, we imposed the constraint that if a field 
overlaps the output image, other good fields from the 
same run must cover the output image RA completely. 

2.2. Photometric Calibration 

The standard SDSS processing calibrates the single 
pass data using data obtained by the 20" PT telescope. 
These include a set of star fields in the Stripe 82 area and 
extinction values measured on the night t he SDSS tele- 
scope data are obtained. The PT pipeline ( | Tucker et al. 
2006 1 calibration results in runs calibrated to rms of 1% 
m gri, 2% i n u and 3% in z, ( Abaz^ianetal.|2003[|Ivezic 



2004|). The DR8 data QAihara et al.|20lTFm cludes 

which 



et al. 

the ubercalibration of Padmanabhan et al.|(|2008[ 



used Apache Wheel scans and runs obtained for Segue. 
The ubercalibration data were not available at the time 
this work was performed. 

Much of the Stripe 82 data was non-photometric so 
we developed a method to calibrate them. In the process 
we also re-calibrated those runs taken under photometric 
conditions. The runs were calib rated following the pre- 
scription of Bramich et al. ( 2008 ) , which builds a catalog- 
level coadd of bright stars to match stars in the frames 
and compute a zero-point shift. The resulting photomet- 
ric calibration of a given run is good to 0.02 mag in up 
to 1 m ag of atmospheric extinction (see also Ivezic et al. 
(2007) for a discussion of calibration through clouds). 

The PT telescope observed the calibration patches in 
Stripe 82 many times while measuring the extinction for 
each standard run independently. Averaging stars cali- 
brated using these independent calibrations will provide 
an increase in photometric accuracy. The increase is 
unlikely to be \A/V as there are systematics floors from 
residual flat field variations and uncorrected atmospheric 
transmission variations. 

We used 62 of the photometric runs for which normal 
SDSS PT calibrations were available to construct a stan- 
dard star catalog. We start with a set of bright, isolated, 
unsaturated stars, with 14 < r < 18, taken from a set of 
high quality photometric runs covering both strips of the 
whole stripe acquired over an interval of less than twelve 
months (2659, 2662, 2738, 2583, 3325, 3388). We then 
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TABLE 1 

Data Selection Criteria 



Scope Criterion 



Description 



Acceptance Rate 



125 < run < 5924 



field r 0.265 v / ncfEpsf < 2.0 
r sky .frames < 375 DN 



r transparency > 1/1.2 

^calibration ^ 1 



data taken on or before 12/1/2005 on Stripe 82 
seeing < 2" 

sky brightness less than 19.5 mag/arcsec 2 

2.5 X the median sky of 150 DN 

allowing at most 0.5 mag increase in sky noise 

less than 0.2 mag of extinction 

enough stars for relative calibration 



91% 
84% 



95% 
95% 



match the individual detections of these stars in each of 
the 62 runs, using a matching radius of 1 arcsec. On 
average, there are 10 independent measurements of each 
star among the 62 runs, and we only include in the ref- 
erence catalog those with 5 or more measurements. We 
then compute the mean of the independent calibrated 
flux measurements of each star and adopt that mean 
flux, defining it separately in each band. We use the 
fluxes measured in the SDSS "aperture 7" , which has a 
radius of 7.43 arcsec; this aperture is typically adopted in 
the SDSS as a reference aperture appropriate for isolated 
bright star photometry. 

Using this standard star catalog, we computed the rel- 
ative zero-point offset of all fields in all runs used in the 
coadd, regardless of whether they were photometric or 
non-photometric runs initially. The relative zero-point 
offset was defined as the median fractional flux difference 
of the standard stars in each field in the run. These field- 
by-field offsets are the atmospheric transmission T, which 
we need for weighting in the coadd as well as to place the 
fields onto the same calibration. There is no requirement 
that T be a smooth function of RA as for example T will 
change with time on non-photometric nights. 

The u band images have signal-to-noise ratio signifi- 
cantly poorer than the other bands and require special 
treatment. All u runs have a provisional calibration ap- 
plied, but in case of non-photometric runs these cali- 
brations are purely the average instrumental zero-point. 
We use this approximate calibration to eliminate lower 
signal-to-noise stars by rejecting those with u > 18. The 
remaining stars are used to match against the standard 
star catalog to find the relative zero-point. 

The flux calibrations are relative magnitude offsets 
from a zero-point, —23.90. We interpret them as vari- 
ations in T with respect to the mean transparency. We 
build a table of linear values of these relative flux scale 
factors T. We T in the coadd image creation process 
to place the images onto the system where the atmo- 
sphere has a uniform transparency, incorporating both 
zero-point and extinction from the all-sky-photometry 
model of the sky. 

TABLE 2 

Runs used in the coadd 
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run MJD date 


RA a tart 


RAcnd 


strip 


phot 


125 51081 1998/09/25 


-10.49 


76.00 


S 


1 


1033 51464 1999/10/13 


-49.00 


-9.40 


N 


1 


1056 51467 1999/10/16 


-35.32 


-0.12 


S 


1 


1752 51818 2000/10/01 


21.45 


79.13 


N 


1 



run 


MJD 


date 


RA s tart 


fiAend 


strip phot 


1755 


51819 


2000/10/02 


-55.68 


47.42 


S 


1 


1894 


51875 


2000/11/27 


31.68 


58.91 


S 


1 


2385 


52075 


2001/06/15 


-53.59 


-37.64 


N 


1 


2570 


52170 


2001/09/18 


17.45 


59.99 


N 


1 


2578 


52171 


2001/09/19 


29.06 


61.44 


N 


1 


2579 


52171 


2001/09/19 


36.46 


60.56 




1 


2583 


52172 


2001/09/20 


-56.42 


-16.95 


s 


1 


2585 


52172 


2001/09/20 


-32.85 


-17.43 


s 


1 


2589 


52173 


2001/09/21 


15.79 


62.58 


N 


1 


2649 


52196 


2001/10/14 


-17.58 


10.29 


N 


1 


2650 


52196 


2001/10/14 


4.18 


31.14 


N 


1 


2659 


52197 


2001/10/15 


-58.26 


-34.58 


N 


1 


2662 


52197 


2001/10/15 


-41.69 


39.74 


TVT 

IN 


1 


2677 


52207 


2001/10/25 


4.25 


39.91 


N 


1 


2700 


52224 


2001/11/11 


20.73 


63.92 


TVT 

IN 


1 


2708 


52225 


2001/11/12 


-15.63 


25.61 


N 


1 


2709 


52225 


2001/11/12 


20.40 


63.25 


S 


1 


2728 


52231 


2001/11/18 


-61.21 


34.24 


N 


1 


2738 


52234 


2001/11/21 


12.86 


62.18 


N 


1 


2768 


52253 


2001/12/10 


-17.40 


35.82 


N 


1 


2820 


52261 


2001/12/18 


20.66 


61.32 


N 


1 


2855 


52282 


2002/01/08 


19.76 


30.66 


N 


1 


2861 


52283 


2002/01/09 


32.84 


66.17 


N 


1 


2873 


52287 


2002/01/13 


14.21 


62.59 


N 


1 


2886 


52288 


2002/01/14 


14.21 


62.52 


S 


1 


3325 


52522 


2002/09/05 


-15.54 


61.16 


s 


1 


3355 


52551 


2002/10/04 


19.37 


61.00 


s 


1 


3360 


52552 


2002/10/05 


-53.45 


25.70 


b 


1 


3362 


52552 


2002/10/05 


20.47 


57.45 


M 


i 
i 


3384 


52557 


2002/10/10 


-54.65 


66.50 


N 


1 


OOOO 


OZOOO 


ZUuZ j 1U/ 11 


A 1 OO 


OZ.OO 


S 


1 


3427 


52576 


2002/10/29 


-51.41 


-24.96 


S 


1 


3430 


52576 


2002/10/29 


20.76 


40.16 


S 


1 


3434 


52577 


2002/10/30 


-52.22 


36.25 


S 


1 


3437 


52578 


2002/10/31 


-50.78 


24.99 


N 


1 


3438 


52578 


2002/10/31 


30.63 


62.61 


S 


1 


3460 


52585 


2002/11/07 


19.41 


61.44 


s 


1 


3461 


52585 


2002/11/07 


42.77 


61.24 


N 


1 


3465 


52586 


2002/11/08 


-34.17 


21.68 


s 


1 


4128 


52908 


2003/09/26 


-7.82 


61.45 


N 


1 


4136 


52909 


2003/09/27 


27.90 


60.79 


S 


1 


4145 


52910 


2003/09/28 


-16.01 


61.85 


s 


1 


4153 


52911 


2003/09/29 


-16.34 


11.87 


N 


1 


4157 


52912 


2003/09/30 


19.23 


61.20 


N 


1 


4184 


52929 


2003/10/17 


-52.92 


-10.08 


N 


1 


4187 


52930 


2003/10/18 


-51.72 


-34.03 


s 


1 


4188 


52930 


2003/10/18 


-15.85 


8.25 


N 


1 


4192 


52931 


2003/10/19 


-52.66 


23.32 


S 


1 


4198 


52934 


2003/10/22 


-53.56 


61.16 


N 


1 


4203 


52935 


2003/10/23 


-60.00 


61.51 


S 


1 


4207 


52936 


2003/10/24 


-55.00 


61.38 


N 


1 


4247 


52959 


2003/11/16 


-15.67 


28.74 


S 


1 


4253 


52962 


2003/11/19 


-15.59 


13.24 


N 


1 


4263 


52963 


2003/11/20 


-16.70 


53.98 


S 


1 


4288 


52971 


2003/11/28 


19.27 


46.65 


s 


1 


4797 


53243 


2004/08/26 


-53.53 


-24.28 


N 


1 


4868 


53286 


2004/10/08 


-30.53 


62.90 


N 


1 


4874 


53288 


2004/10/10 


-62.40 


88.08 


N 


1 
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TABLE 2 
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run 


MJD 


date 


R A * * 


R A j 


strip phot 


4895 


53294 


2004/10/16 


-4.80 


70.77 


N 


1 


4905 


53298 


2004/10/20 


0.38 


72.29 


N 


1 


4917 


53302 


2004/10/24 


-65.59 


52.81 


N 





4930 


53313 


2004/11/04 


-58.01 


1.81 


S 


1 


4933 


53314 


2004/11/05 


-53.70 


63.36 


N 


1 


4948 


53319 


2004/11/10 


7.31 


62.39 


N 


1 


5042 


53351 


2004/12/12 


18.42 


61.43 


S 


1 


5052 


53352 


2004/12/13 


-15.67 


25.72 


S 


1 


5566 


53616 


2005/09/03 


-33.61 


60.28 


N 





5582 


53622 


2005/09/09 


-55.62 


58.95 


S 





5590 


53623 


2005/09/10 


-60.69 


12.78 


N 





5597 


53625 


2005/09/12 


-64.68 


-17.02 


S 





5603 


53626 


2005/09/13 


-66.47 


63.02 


N 





5607 


53627 


2005/09/14 


-63.90 


62.25 


S 





5610 


53628 


2005/09/15 


-66.73 


65.38 


N 





5619 


53634 


2005/09/21 


-64.44 


63.20 


S 





5622 


53635 


2005/09/22 


-64.48 


63.38 


N 





5628 


53636 


2005/09/23 


-64.61 


21.61 


S 





5633 


53637 


2005/09/24 


-61.59 


59.71 


N 





5637 


53638 


2005/09/25 


-21.61 


63.42 


S 





5642 


53639 


2005/09/26 


-10.71 


62.83 


N 





5646 


53640 


2005/09/27 


-65.64 


70.65 


S 





5658 


53641 


2005/09/28 


15.19 


56.03 


N 





5666 


53643 


2005/09/30 


40.36 


63.58 


S 





5675 


53645 


2005/10/02 


-60.60 


-40.41 


S 





5681 


53646 


2005/10/03 


22.34 


54.21 


s 





5709 


53654 


2005/10/11 


-67.46 


24.71 


N 





5713 


53655 


2005/10/12 


-68.40 


42.80 


s 





5731 


53657 


2005/10/14 


20.28 


62.34 


N 





5732 


53657 


2005/10/14 


46.22 


62.33 


S 





5754 


53664 


2005/10/21 


-57.68 


59.33 


S 





5759 


53665 


2005/10/22 


-59.55 


59.24 


N 





5763 


53666 


2005/10/23 


-59.22 


5.29 


S 





5765 


53666 


2005/10/23 


1.41 


56.19 


N 





5770 


53668 


2005/10/25 


-56.57 


59.20 


N 





5771 


53668 


2005/10/25 


32.20 


62.16 


S 





5776 


53669 


2005/10/26 


-59.98 


59.26 


S 





5777 


53669 


2005/10/26 


23.52 


59.29 


N 





5781 


53670 


2005/10/27 


-56.20 


59.17 


N 





5782 


53670 


2005/10/27 


31.77 


62.19 


S 





5786 


53671 


2005/10/28 


-36.66 


63.30 


S 





5792 


53673 


2005/10/30 


-62.44 


59.30 


N 





5797 


53674 


2005/10/31 


-59.00 


59.39 


S 





5800 


53675 


2005/11/01 


-59.49 


59.98 


N 





5807 


53676 


2005/11/02 


-48.55 


59.20 


S 





5813 


53677 


2005/11/03 


-65.13 


45.97 


N 





5820 


53679 


2005/11/05 


-45.43 


62.22 


S 





5823 


53680 


2005/11/06 


-60.09 


62.32 


N 





5836 


53681 


2005/11/07 


-60.56 


62.35 


S 





5842 


53683 


2005/11/09 


-59.31 


62.22 


N 





5847 


53684 


2005/11/10 


-63.50 


63.90 


S 





5866 


53686 


2005/11/12 


17.32 


63.41 


N 





5878 


53693 


2005/11/19 


-62.67 


64.29 


N 





5882 


53694 


2005/11/20 


-63.69 


63.13 


S 





5889 


53696 


2005/11/22 


35.34 


63.13 


S 





5895 


53697 


2005/11/23 


-63.25 


62.62 


S 





5898 


53698 


2005/11/24 


-65.81 


62.31 


N 





5902 


53699 


2005/11/25 


-62.81 


62.75 


N 





5905 


53700 


2005/11/26 


-68.06 


62.74 


S 





5918 


53704 


2005/11/30 


-62.39 


62.28 


N 





5924 


53705 


2005/12/01 


-63.20 


62.51 


S 






3. COADD IMAGE CREATION 

We aim at coadding the Stripe 82 data and running 
it through PHOTO. To make the coadd we need the data 
images and weight maps. We build a weight map by mul- 
tiplying an inverse variance map and a geometry mask. 
PHOTO requires a map of the saturated pixels. In this 
section detail the process of creation of each of these im- 
age components and describe how we use them in the 
coaddition. 



3.1. Sky Subtraction 

The sky brightness varies both spatially and tempo- 
rally, night to night and due to clouds. We removed the 
sky before mapping. As the data we used for the coadd 
came from DR7 and earlier, the improved sky subt raction 
of the curre nt DR8 PHO TO was not implem ented (Aihara 
et al. (20111; see also Blanton et al. (2011) for continued 

ork on this subtle problem). The DK7 PHOTO algorithm, 
which was current at the time, produced a sky image by 
calculating the median of the 256 x 256 pixel boxes in 
the data image on a grid of 128 pixels in each dimension. 
The sky was then determined using bilinear interpola- 
tion. This algorithm over-subtracts the extended parts 
of galaxies on the scale of the 128 x 128 pixel grid used in 
the sky calculation, leaving artifacts due to astrophysical 
objects in the data. PHOTO sky subtraction engine was 
therefore deemed not suitable for our purposes and we 
developed our own method. 

Subtracting a global sky for each frame would work 
poorly as the sky changed with time and thus with row 
number in the SDSS data. Instead we adopted a sky 
value that was allowed to vary linearly with row number 
and thus time. For a given frame and the two frames on 
either side of it in the run (and after removing the SDSS 
soft bias of 1000 DN), the median of the 2048 pixels 
along each row was calculated, resulting in a 3 x 1489 
pixel sky vector. This vector is a time series estimate. 
In order to deal with bright stars, the rms of the vector 
was calculated using a 3er 5 iteration sigma clipping, and 
the rms was used to reject pixels more than 2a from the 
mean. A linear least squares fit to the remaining vector 
was used to model the sky, which was subtracted from 
the image row by row. 

3.2. Astrometry 
The coadd relies on the existing astromet ry produced 



using the astrom pipeline (Pier et al. 2003). All of the 



images used in t he coadd had astrometric calibrations 



Pier et al. ( 2003 1 were able to achieve positions accurate 
to ~ 45 mas rms per coordinate by calibrating to the US 
Naval Observatory CC D Astrograph Catalogue (UCAC; 
Zacharias et al. ( 2000[ )). The accuracy is limited primar- 
iiy by the accuracy of the UCAC positions (~ 70 mas 
rms at the UCAC survey limit of R « 16) and the density 
of UCAC sources. This accuracy can be represented in 
the affine transformations that are standard in the WCS 
convention. 



Pier et al. 



(2003) also noted that there are 
systematic optical distortions due to the camera present 
in the data. We will use the astrom measurements to 
remove these distortions. 

The process of forward mapping requires a transfor- 
mation from RA,Dec (a,S) to pixel location in the input 
image. The SDSS runs were taken along great circles. 
Thus astrom worked in a coordinate system in which 
each run's great circle is the equator of the coordinate 
system. In this great circle coordinate system, the lati- 
tude of an observed star never exceeds about 1.3°; thus 
the small angle approximation may be used and lines of 
constant longitude are, to an excellent approximation, 
perpendicular to lines of constant latitude. Longitude 
and latitude in great circle coordinates are referred to as 
fi and v, respectively, v is equal to along the great cir- 
cle, y, increases in the scan direction, and the origin of \x 



G 
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is chosen so that [i = 02000 at the ascending node (where 
the great circle crosses the J2000 celestial equator). The 
conversion from great circle coordinates to J2000 celestial 
coordinates is then 

sin(/i — Mo) cos v cos i — sin v sin i 

tan(a-/i ) = -, r (1) 

cos(/i — Ho) cos v 

sin {$) = sin(/i — fi ) cos v sinz + sin v cos i (2) 

where i and fio are the inclination and J2000 right as- 
cension of the great circle ascending node, respectively. 
fin = 95° for all survey stripes, and for Stripe 82 i w 0. 

Given the great circle coordinates (jx, v) we can trans- 
form to distortion corrected frame coordinates {x 1 ,y') us- 
ing the affinc transformation 



fJ-CMP = a + bx' + cy' 
Vcmp = d + ex' + fy' 



(3) 
(4) 



The transformation from (a/, y') to {x, y) accounts for 
optical distortions which, in drift-scan mode, are a func- 
tion of column only: 



V = V 



go 

h 



9iy + 92y - 

h\y + h 2 y 2 



93V 
-h 3 y 3 



(5) 
(6) 



Equations fl]6]) provide the pixel coordinates on the in- 
p ut image that correspo nds to a given (RA,Dec) position. 



Bramich et al. ( 2008J) performed a recalibration of the 
Stripe 82 astrometry using a run from the mid-time of the 
Stripe 82 observations as a reference. This removed an 
erroneous but measurable galaxy mean proper motion of 
~ 10 mas/yr in both RA and Dec due the p roper motion 
of reference stars (see Bramich et al. ( |2008[ ), Fig. 4). The 
optics distortions removed in this sec tion have a maxi - 
mum peak to peak shift of 80 mas (see Pier et al. (2003), 
Fig. 2) and are larger than the astrometry shifts induced 
by stellar motion over the 5 year range of the data used 
in the coadd. Therefore, although ideally we should have 
removed the reference star proper motion drift, in prac- 
tice these have little effect. 

3.3. Astro-metric Mapping of Data Images 

We geometrically map the input images onto the out- 
put image. We defined output frames aligned along the 
J2000 equator with rows aligned perpendicular to RA, 
in a standard SDSS image format, from —50° < a < 60° 
and -1.25° < S < 1.25°. 

Since the output image is simply a locally flat tangent 
projection of the sky, the mapping must remove optical 
distortions and provide a surface brightness estimate at 
the aligned pixel locatio n. To perform the mapping we 
used a version of Swarp (Bertin et al.l 2002 1 modified to 



perform the astrometric conversions described in {3.2 



Each pixel in the mapped image is estimated from the 
input image pixels using a Lanczos interpolation kernel, 
which is a truncated sine interpolation. Given bandwidth 
limited signal of infinite extent, sine function interpola- 
tions reproduce exactly the data after resampling. Our 
data is not undersampled (it has 0.4" pixels and seeing 
of ~ 1-3"), but it is not of infinite bandwidth either and 
this motivates the use of a truncated kernel. 

We used a two dimensional Lanczos-3 kernel retain- 
ing 3 maxima on each side of the center in each dimen- 
sion. The one dimensional Lanczos-3 kernel is L(x) = 



sinc(x)sinc(x/3), for —3 < x < 3. 
sional interpolation formula is 



Then the 2 dimen- 



I(r,c)=J2Hi,j)L(r-i)L(c 



J) 



(7) 



1,3 



where r, c are the output image pixel coordinate and i, j 
are the input image pixel coordinate. The Lanczos-3 win- 
dow is well-behaved in terms of reduction of aliasing, 
minimal ringing, and lack of smoothing, but a Lanczos-3 
interpolation does use (2a;3 + l) 2 = 49 pixels to estimate 
the value of one output pixel. This is too large to use 
for bad pixels (e.g, saturated pixels), so we use nearest 
neighbor interpolation and reduce the weight of bad pix- 
els duri ng L anczos-3 interpolation by using a mask (see 
section { 3.8 1. 



3.4. Inverse Variance Map 

To keep track of the variance of data images, pixel by 
pixel inverse variance images are a natural choice, but 
they produce biases in the resulting mean. At low sig- 
nal to noise the upward fluctuations in signal are given 
more weight than downward fluctuations as a result of 
the one-sided nature of the Poisson distribution. This 
bias is deterministic and one could correct for it, but, 
for example, for u-band data with its 120e~ of sky noise, 
pixels at la above sky would be biased by 0.5% and this 
would be a fair fraction of our photometric error bud- 
get. Another problem is that per-pixel inverse variance 
weighting systematically changes the shape of the PSF 
as a function of the magnitude of the object. This would 
cause serious complications to our PSF-based photome- 
try using PHOTO. For these reasons we chose a different 
method. 

We computed the variance of the sky as measured on 
the frame from the width of the sky histogram. As we 
used a linear gradient sky subtraction, each image is as- 
signed a variance image that is the variance of the sub- 
tracted linear sky gradient. This of course assumes Pois- 
son statistics while the data are in ADU. In calculating 
this variance we did not include the effective gain, g e g, so 
gain variations are not accounted for. These variations 
are < 30% though. 

3.5. Geometry Mask 

The geometry mask, keeps track of which pixels in the 
input image actually contribute to the coadd. In this 
mask definition we account for image defects found by 
PHOTO, in particular, cosmic rays, saturated pixels and 
bad columns. 

PHOTO produces a 16 bit mask image, the f pM file. For 
the geometry mask we are interested in the INTERP bit. 
The INTERP flag is set for any pixel for which PHOTO 
used interpolation to fill its value. This happens for cos- 
mic rays, saturated pixels, and bad columns. These pix- 
els are poor but not useless estimates of the true value 
of the pixel. We set the geometric weight of such pixels 
to a small value, 0.0000001. This ensures that if there is 
no input image which contributes a good estimate for a 
given output frame pixel, such as the center of a satu- 
rated star, the interpolated value is used. 

The geometry of the SDSS images are also encoded 
into the geometry mask. In SDSS images, the first 124 
rows of each frame are duplicates of the last 124 rows of 
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the preceding frame. This replication of pixels was done 
so that objects could be well measured despite being on 
edge of a frame. To account for this, the first 124 rows of 
the mask are set to 0. In the SDSS images there are also 
scanline to scanline overlaps between North and South 
strips. This is actual exposure time, again designed to 
make objects on the edges measurable. For our purposes 
we kept the North and South strip data separate as the 
extra exposure was on too thin a strip to be useful. 

3.6. Satur Mask 

PHOTO's bit mask also contains information about sat- 
urated pixels, encoded in the SATUR bit. The SATUR 
flag is set for any pixel that PHOTO determines to be sat- 
urated. We set these pixels to 1 in an image otherwise 
filled with zeros. This is the SATUR map, which we need 
for running PHOTO on the coadded images. 

3.7. Weight Map 

Although closely related, the weigh map and the in- 
verse variance image are not the same. The weight map 
is the inverse variance image multiplied by the geometry 
mask. We set all masked non-zero pixels to the INTERP 
flag value (0.0000001). This allows us to keep track of 
the pixels altered by the masking. 

We use the weigh map as input for the coaddition pro- 
cess, in which a weighted clipped mean of the data im- 
ages, all mapped onto the same output image, will be 
performed. In addition, we coadd the inverse variance 
images and the satur maps as well, and although for those 
we use straight sums, the weigh map plays a role in the 
sense that only the pixels that pass the clipping for th e 
data coadd are included (for details, see section £ 3.9.2). 



3.8. Astrometric Mapping of Map & Mask Images 

The inverse variance map, satur mask and weight map 
are all mapped onto the same output image as the data. 
For the inverse variance and weight images we apply the 
same Lanzcos-3 interpolation used for the data in or- 
der to replicate the noise correlation. The satur masks, 
propagating saturated pixels, were mapped using a near- 
est neighbor interpolation. This is more suitable than 
Lanzcos-3 because masked pixels an area of influence lim- 
ited to their first neighbors. 

In the process of performing the wide Lanzcos-3 in- 
terpolation of the data, the geometry mask is used to 
prevent masked pixels from contributing to the interpo- 
lation with more than minimal weight. In addition, once 
the images are mapped onto the output image, any non- 
overlap is set to zero weight using the geometry mask. 

After this mapping procedure, each stack of images is 
ready for the coaddition process. No scaling has been 
done at this stage, so the stack for a given output image 
can be further filtered and/or weighted as needed. 

3.9. Coaddition 

With the four stacks of images (data, inverse variance, 
weight and satur) all aligned and cropped to match the 
output image, we have all the inputs needed to produce 
the coadded images and run PHOTO on them. We chose 
to use all the data in our image set. This allows us to 
maximize the depth of the coadd data. 

Alternatively, one could devise a selection criteria to 
produce a coadd dataset tailored for a specific purpose. 



For example, one could select only images with the best 
seeing and limit them in number to form a uniform depth, 
aiming at weak lensing studies. 

3.9.1. Weights 

The data in our image set is of variable quality and we 
wish to optimize our coadd, so we designed a weighting 
scheme. The signal to noise ratio of the measurement of 
flux from a star is: 



S/N oc 



N, 



photons 



V ^psf^sky 



(8) 



where iVphotons is the number of photons detected from 
the source, A ps f is the area in pixels the source subtends, 
and <7 s k y is the sky noise per pixel. As N is propor- 
tional to transparency T and A cx FWHM 2 , we use the 
following as weights: 



Ti 



FWHMjV^ft) 



(9) 



This gives highest weight to good seeing data taken when 
the sky is clear and dark. We choose these weights be- 
cause although the usual inverse variance weighting pro- 
duces the minimum variance image, the signal to noise 
for a star depends on square of seeing. With this weight, 
the PSF of the coadd is 0.3" less than the median of the 
input images in ugri and 0.2" less in z. 

We take the seeing for each field and filter from the 
tsField file, which is simply an average over the frame, 
and use it regardless of how much the frame contributes 
to a given output coadd frame. The transparency T, 
in turn, is already available as a product of the pho- 
tometric calibration discussed in section 92.21 We have 



therefore cataloged for each image the transpareancy, the 
seeing, and via the inverse variance image, the sky vari- 
ance. This allows us to proceed to the coaddition. 

3.9.2. Weighted Clipped Mean 

We coadd by building the stack of mapped images on 
a given output frame and performing a weighted clipped 
mean. The weight images described in |3.7| are multiplied 
by T/FWHM 2 . Flux calibration of our data images by 
1/T means that our weight images get another factor of 
T 2 due to error propagation. 

Then for each output pixel we collect the corresponding 
mapped pixels, and reject outliers, using a iterative 5cr 
rejection, where for "c" we use the 25-75% interquartile 
range The final average uses the weights of Eq. [9| 

The data images are coadded as described above. The 
inverse variance image and the satur mask coadds were 
computed using straight sums (not averages) , though us- 
ing only those pixels corresponding to those that make it 
past the clipping process into the coadded data image. 

3.10. Output Images 

Figure [2] shows a side by side comparison between 
coadd single pass data in r-band. The single pass coun- 
terpart is one out of 28 images used in the coaddition 
which means that it passed all the cuts discussed in the 
text. This example features run 206, camcol 3, field 505 
in r-band. It illustrates the fact that a number of ob- 
jects bellow the detection threshold of each image can be 
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well detected and measured in the coadd. The seeing, 
however, is larger for the coadd image in this example, 
indicating that this particular single pass image has see- 
ing better then the median seeing on the stack. 

4. CATALOG CREATION 

The coadd production yields 800 fields along six 
columns in each of two stripes and each field includes 
data in five bands. The resulting 48,000 images are pro- 
cessed through a modified version of the SDSS image 
processing pipeline PHOTO, to yield object catalogs. 

4.1. PSF Measurement 

Our principal challenge is to measure the PSF of the 
coadd. The PSF of SDSS images varied in time, corre- 
sponding to image rows, due to atmospheric fluctuations. 
It also varied in space, corresponding to columns, due 
to camera optics. In the standard single run reduction, 
PHOTO fits the spatial variation of the PSF by computing 
a PSF basis set using a Karhunen-Loeve (KL) expan- 
sion. The stars in the frame and the two flanking on 
either side (five in total) are used to determine the KL 
basis functions Br(u,v): 



Pi(u,v) 



,a r iB R {u,v) 



(10) 



where Pi is the PSF of the i th star, u, v are the pixel 
coordinates relative to the basis function origin, and n 
sets the number of terms to use in the expansion (we use 
n = 3). The stars in the image and the flanking half an 
image on either side (two in total) are used to determine 
the KL coefficents aF: 



l+m<N 
l=m=0 



(11) 



where x, y are the coordinates of the center of the i th 
star, N is the highest power of x or y included in the 
expansion (we use N — 2, a quadratic spatial variation 
of the PSF), and the b r lm are found by minimising 



(Pi(u,v) 



alB r (u,v)f 



(12) 



In the coadd we took advantage of the fact that the KL 
basis set and coefficients had already been determined 
each input image as each of the input runs had already 
been separately processed through PHOTO. Using the KL 
basis set for each input run corresponding to the given 
output frame and the two flanking frames, we compute 
a weighted sum of the model PSFs (using the weights 
of equation M on a 2-dimensional grid with a spacing 
of ~ 1.5'. We then fit the KL basis -Br(u,i;)'s to these 
over the coadd frame and the two flaking frames. The 
computation of the coefficents of the PSF expansion a\ 
are done on the coadd frame and the two flanking half- 
frames as in the single pass runs, with the exception that 
we set the maximum number of coefficients, n in Eq. |10[ 
to n — 4 and highest power of x, y, N in Eq. [TT] to 
N = 3. 



4.2. Effective Gain & Sky 

PHOTO had to be modified to read in a file containing 
the weights and effective gains of the images. The ef- 
fective gain of a coadd image G should be such that it 
gives bright objects Poisson statistics for the variance as- 
sociated with a source when scaling the averaged pixel 
counts back to electrons. Therefore, 



a 2 (P) = P/G 



(13) 



where P is the weighted sum, over all exposures i, of the 
object counts Pi, 



p=zl Wipi = E wiTi0i 



(14) 



and pi = TiOi is the object's flux per exposure Oi times 
the flux scale factor Tj. The corresponding variance is: 



a 2 (P)- 



■ ^2wfa 2 (pi) 
^wfTfc^ioi) 
■^wfTfoi/gi 
■y^PiWiTi/gi 



(15) 



Assuming that pi does not vary much from image to im- 
age, we obtain: 



\P) « (ft) £ wlTi/gi = P^w^/gt 



(16) 



Eqs. [13| and [l6| we conclude that the effective gain is: 



G= [j2^i T i/9i 



(17) 



In the SDSS the gains gi of all frames going into a single 
coadd frame are the same (it is always the same CCD), 
so they may be replaced by a single g outside the sum: 

9 



G 



(18) 



To process the images with PHOTO we also need to 
compute the effective sky level S' . Assuming that dark 
and read noise are accounted for in the Poisson noise 
of the effective sky, we can easily obtain S' using the 
effective gain G calculated in Eq. [18] and the relation 
S" = Ga 2 (S'), which is analogous to Eq. 



13 



We used 

the clipped variance of the coadd frame to compute 
er 2 (S"). Alternatively, one could take the sky variances 
and weights of the input frames and compute it using 
v 2 (S') = Ew 2 T 2 a 2 (s t ). 

Using this method we obtained, for each field and for 
all 5 filters, the effective sky (S"), sky noise (a(S')) and 
gain (G, Eq. 18). With these quantities in addition to 
the weights (wi, Eq. [9]) we processed all of the coadd 
fields through PHOTO. 

4.3. Applying the Calibration 
The SDSS code Target is used to apply the calibration 



to the raw outputs of PHOTO. Following |Lupton et al. 
(19991, we convert from fluxes / to mags m using 



m(f/f ) 



2.5 
In 10 



asinh 



f/fo 
2b 



■ln(6) 



(19) 



where b is an arbitrary dimensionless softening parameter 
below which the magnitude scale goes from linear to loga- 
rithmic and /o is a reference flux that sets the zero-point, 
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Fig. 2. — Comparison between single pass (left) and coadd (right) images in r-band for run 206, camcol 3, field 505, RA=15, Dec=0. 
Images are shown with the same scale, contrast and stretch. The single pass counterpart (run 5800, camcol 3, field 505) is one out of 28 
images used in the coaddition of this particular image. This example illustrates the fact that a large number of objects bellow the detection 
threshold of each image can be well detected and measured in the coadd. 

exponential profile convolved with the PSF) and r ps f is 
the PSF magnitude at the position of the object. This 
simple estimator performs well for single pass data: 95% 
correct at magnitude 21 in r. 

In the case of the coadd data, star count plots showed 
dramatic increases in the numbers of stars at magnitudes 
where galactic models do not. This suggested that more 
numerous galaxies were being misclassified as stars there. 
Thus we instead used the more stringent criterion 



TABLE 3 
Asinh Magnitude Softening 
Parameters for the Coadd 



filter 


b 




zp 


ro(10&) 


u 


1.0 x 10- 


ii 


27.50 


24.99 


9 


0.43 x 10" 


-11 


28.42 


25.91 


r 


0.81 x 10" 


-11 


27.72 


25.22 


i 


1.4 x 10- 


-11 


27.13 


24.62 


z 


3.7 x lO- 


-11 


26.08 


23.57 



Note. — Values reported by |Abazajian| 
|et al.M2009) . Column zp is the zero-point 
magnitude, zp = m(0). The final column 
gives magnitude associated with an object 
for which f/fo = 106. 



zp = m(0), of the magnitude scale. The values of b that 
we used for the coaddition is given in Table [3j along with 
the asinh magnitudes associated with a zero flux object 
and the magnitudes corresponding to / = 10/o&. Above 
this scale, the asinh magnitude and the traditional loga- 
rithmic magnitude differ by less than 1% in flux. These 
values can be compared to their equivalent numbers for 
the m ain survey, given in Table 21 of Stoughton et al. 
(2002). The coadd images were all placed onto a uni- 



Fpsf 



□dell < 0.03 



(20) 



to select stars. As PHOTO measures every parameter for 
every object regardless of its determination of object type 
this has no effect on the other measurements. 

5. DATA PRODUCTS VERIFICATION 
5.1. Photometric Calibration 



We use the standard star catalog of Ivezic et al. (20071 
to v erify our photometri c calibration. To build that cat 



form flux scale such that 1 DN corresponds to a flux 
of 1 picomaggie, corresponding to a logrithmic (and not 
asinh) mag of 30. 



4.4. 



Star/Galaxy Separation 

The SDSS standard star/galaxy separati on simply 
classifies as stars all objects in the region (Abazajian 
et al.||2004"|) |r psf - r mode i| < 0.145 where r m dei is the 



model magnitude (the best fit galaxy deVaucouleurs or 



alog Ivezic et al. (2007) took the median of individual 
measurements of bright stars from 58 Stripe 82 runs and 
then applied several corrections to their catalog: 1) a 
color-term like correction for the bandpass of each cam- 
era column; 2) a purely RA flat field correction in r de- 
rived by comparing PT data with SDSS data; and 3) a 
purely Dec flat field correction to the colors relative to 
r band from stellar locus colors. We will take the re- 
sulting catalog (which is calibrated to 1% accuracy) as 
the truth and com pare with our own measurements. As 
Ivezic et al. (2007) applied corrections that we did not, 
the comparison is not completely circular. 

We sliced our star cat alog into a series o f 1 magnitude 
bins and matched to the Ivezic et al. ( 2007 ) catalog using 
a 1" matching radius and discarding objects with more 
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than one match inside that radius. No flags were applied 
to our star selection; in particular we did not demand an 
isolated, well measured set of stars to start with. 

Table [4] summarizes the results of our comparison. We 
defined Aj as the median of the difference between our 
measurements and the quantities reported in the stan- 
dard star catalog, for magnitudes (i = m) and colors 
(i = c). Ai is a measure of the zero-point offset. Statisti- 
cal uncertainty in the zero-points, obtained as the rms of 
the differences, are of the order a few millimag/ y/N hj ; as 
usual in photometry we can expect systematics to dom- 
inate this. The offsets from the standard zero-points are 
less than 5 millimags (1% photometry corresponds to 10 
millimags) in all cases. 

We also examined the spatial variations of the zero- 
point offset and its uncertainty. To this end, we first 
took 50 equally spaced bins of width 2.2° in RA, and 
computed the mean differences in magnitudes and colors 
for each bin. We repeated this procedure in Dec bins, 
choosing 30 bins of 5' width. Figure [3] shows these mean 
zero-point offsets as a function of RA and Dec, indicat- 
ing that spatial variations non-negligible, specially for u 
band. As a measure of the overall offsets we computed 
the median of those means, (Aj)^, where j means either 
RA or Dec. These values are also included in Table|4]with 
the uncertainties estimated as the rms or the means. 

Our comparison to the standard star catalog reveals, 
therefore, that along the RA axis in Stripe 82 the cali- 
bration varies 5 millimag in g,r,i,z, g — r, r—i, i — z< 10 
millimag in u and u— g. Likewise, along the Dec axis the 
calibration varies by < 5 millimag in g,r,i,z, g — r, r — i, 
i — zand 30 millimag for u, 20 millimag for u— gand < 10 
millimag for the others. These can be seen in Fig. [3] (see 
also Table [IJ). The variation in Dec is significantly larger 
than that m RA, probably reflecting systematic errors 
in the PT flat field images and thus in the PT standard 
overlap fields the used for the calibration. 

5.2. PSF Modeling 

We verified how well we modeled the PSF in several 
ways, starting by computing the spatial variations in see- 
ing, which can be written as 



FWHM(arcsec) = 2^2 ln2 v / mrrcc/2 ■ S (21) 

where S — 0.396" is the pixel scale and mrrcc is the sum 
of the second moments of the PSF. We used high S/N 
stars in each bandpass for this test. Our results, illus- 
trated on the left panel of Fig. |4j show that the seeing is 
best in the redder filters, consistently with a Kolmogorov 
seeing law with the exception of r-band and i-band. We 
interpret this as an effect of selecting the input frames 
for the coadd in r-band (see Table [T]) associated with an 
effect of the time scales of the Kolmogorov law. Recall 
that the data making up a field are taken at different 
times, ranging over 8 minutes from r-band to g-band. 
Our data support the assumption that the time scales of 
Kolmogorov seeing are such that the seeing is uncorre- 
lated after ~ 1 minute, but since i and r are next to each 
other in the imager, they can be correlated. The median 
of the seeing of the single pass images are a few tenths 
of an arcsecond worse than the points on this plot; this 
reflects the weighting of the coadd (Eq. [9]) . 

Figure [4] also show the seeing as a function of Dec, 
averaging over RA. The seeing is affected by the camera 



optics, which causes the upturn at one end of the camera, 
The N strip (run 206) has about 0.075" worse seeing in 
ugiz, and 0.15" worse seeing in r, then the S strip (run 
106), presumably due to the statistics of the seeing in the 
input images. However, at Dec > +0.5°, as the camera 
optics begin to dominate, the seeing difference becomes 
negligible. 

Another interesting test of our PSF modeling is to 
check the reconstructed PSF at the position of stars. In 
Fig. [5]we show the mean ratio of mrrcc for stars and the 
reconstructed PSF at the locations of stars, as a function 
of Dec. The strong declination dependence seen in Fig. [2] 
is not apparent. This suggests that the modeling is fit- 
ting the spatial variation of the PSF well. We also found 
no correlation between the statistics of star galaxy sepa- 
ration and column number, again suggesting reasonable 
success in the PSF modeling. 

However, given the importance of the PSF modeling 
for crucial aspects of the data, such as star/galaxy classi- 
fication, accurate photometry and shape measurements, 
we examine it in more detail. A sensitive test of the 
PSF modeling is the r-band "PSF minus model" plot, 
r psf — '"model vs - r psf ■ This is shown in Fig.[6]for both stars 
(red) and galaxi es ( blue) (for details on the star/galaxy 
separation, see [4.4|. Stars are expected to be found at 
a very narrow region around r ps f — r mo d G i = 0. and this 
is clearly the case for the bright magnitudes in our plot. 
There is a noticeable trend towards negative values of 
fpsf _ f*modoi = at the faint end, around magnitude 
23, and an upward spread at about magnitude 22. The 
upward spread cannot be due to galaxies behind stars 
causing the PSF to slightly broaden as this was excluded 
by checking the spatial statistics. These two features in 
the diagram suggest that we have magnitude dependent 
PSF fitting problems which may introduce systematics 
in analyses involving the coadd data. 

While probing this, we found that during the process- 
ing described in section 



the error propogation va 



3.9.2| the coadd incorrectly set 



ue to T instead of the correct 



value of T 2 . This means that the weight actually used 



was: 



1 



w 



FWHMV 2 



(22) 



This has little effect on the coadded image as we only 
selected those images with l/T < 1.2. The PSF pho- 
tometry, however, depends so sensitively on the modeled 
PSF that this may be the root of the problem. The 
PSFs wer e constructed using the correct weights (sec- 
tion ^4.1), but these differed from those of the images. 



In an experiment we ran Photo with input weights that 
were pure inverse variance instead of SN weight that the 
images were built with. The resulting psf-model vs psf 
plots showed stars with deviations away from zero at 0.05 
mags/mag level. The effects we are seeing in Fig. p^are 
at the ~ O.Olmags/mag level, consistent with the lower 
affect the tranmission has on the weights. We conclude 
that the low level psf problems that the structure in Fig. [6] 
indicates is due to this mismatch of weights between the 
coadd images and the coadd photometry measurements. 

5.3. Galaxy Catalog Purity 

These misclassified stars are an issue as they lie in 
the magnitude-color space of LRGs. The contamination 
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TABLE 4 

Relative Zero-points for Selected Coadd Catalog Samples 



filter 




mag 


A m 


color 


A c 


(A m > R A 


(Ac)ra 


(A m )Dec 


(A c )Dec 


u 


1311 


20-21 


1.2 ± 1.6 


u-g 


-4.3 ±3.6 


-21.0 ± 8.2 


-12.1 ±8.1 


-21.3 ± 28.8 


-12.4 ±15.3 


g 


1399 


19-20 


2.2 ± 1.6 


8-1' 


-1.7 ± 1.4 


6.5 ±4.6 


-1.2 ±3.5 


5.8 ±10.2 


-1.1 ± 6.3 


r 


3703 


19-20 


3.7 ± 1.5 


r-i 


-2.9 ±2.2 


3.2 ±2.5 


4.0 ±2.1 


3.2 ± 6.3 


4.0 ± 4.1 


i 


7436 


19-20 


0.9 ±2.0 


i-z 


-1.5± 2.3 


0.1 ± 2.0 


6.3 ±3.1 


1.0 ± 6.2 


5.7 ± 6.8 


z 


3963 


18-19 


4.3 ±2.4 






-5.7 ± 3.3 




-5.3 ± 11.7 





Note. — Rows correspond to samples created independently for each filter. Magnitude ranges are indicated by 
the mag column. N ^j is the number of stars in each sample. A; is the median zero-point difference, in millimags, 
for either magnitudes (i = m) or colors (i = c) and (Aj) j is the median of the mean difference in spatial bins, RA 
or Dec. 
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Fig. 3. — Zero-point offsets measured in bins of RA and Dec, showing spatial variations in the photometric calibration of the coadd 
data. 



level of the galaxy catalog by missclassified stars was esti- 
mated using the morphological-independent redshift sur- 
vey catalog from the V irmos VLT Deep Survey (VVDS, 



Le Fevre et 
2^hr field 



al. 



2005 1, in particular the the i < 22.5 
We selected the objects above 95% confidence 
level, which resulted in a catalog containing 5158 stars 
and 4264 galaxies. This catalog served as a truth table, 
to match galaxies in the coadd at RA ~ —25°. 908 spec- 
troscopic stars and 3438 spectroscopic galaxies matched 
galaxies in the coadd catalog, indicating a contamination 
level of 18%. This is ~ 3 times higher than expected, as- 
suming that the performance of the SDSS Star/Galaxy 
separator would be as good for the coadd as for the sin- 
gle pass data. The field in question is among the low- 
est galactic latitudes and highest stellar densities in the 
coadd, so this contamination rate is approximately an 
upper limit. The actual contamination rate would be 



roughly proportional to RA in the coadd. 

We verified that most of the problematic objects are in 
a localized region of the psf— model vs. magnitude space, 
specifically inside the triangular region 

Vf - r modcl < 0.1 • r psf - 2.1 

~ r modcl > 0.03 (23) 
21 < r psf < 22.5 

This indicates that an improved star galaxy separator 
can be designed using morphology and presumably color 
cuts. 

5.4. Star & Galaxy Catalog Completeness 



We used the package 2DPHOT ( |LaBarbera et al.|2008 1 
to measure the completeness of the coadd star and galaxy 
catalogs. This is done by adding simulated objects to the 
images and computing the recovery rate for both stars 
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Fig. 4. — Left: Mean FWHM as a function of RA for the five filters. The range of seeing values is consistent with the expected 
Kolomgorov A -0 ' 2 scaling with the exception of r and i bands. As the input images were selected in r and data in the other bands are 
taken at a few minutes removed in time in the order riugz, one might expect this b ehavior. Right: Mean F WHM as a function of Dec. At 
Dec > 0.5 the seeing gets worse as the camera optics begins to dominate (see, e.g., Stoughton ct al. (2002|). 

and galaxies. The parameters used to generate the ob- 
jects are taken from the image itself. 2DPHOT first de- 
tects the objects in the image, performs star-galaxy sepa- 
ration and measures the photometric and structural (Ser- 
sic) parameters (PHOTO is not used in this process). Then 
it creates a list of objects that reproduces the magnitude 
and size distributions of the stars and galaxies found in 
the image and adds these simulated objects to the im- 
age. Finally, it measures the new image and computes 
the completeness as the fraction of objects recovered in 
each magnitude bin. 

The resulting 2DPHOT completeness vs. magnitude 
curves, C(m), are well fit by a Fermi-Dirac distribution 
function 

C ( m ) = i , Z^A ..w^ ( 24 ) 
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- The declination projection of the ratio of size, mea- 
sured through m r r c c, of stars and the psf evaluated at the star 
positions (mrrcc/mrrcc_psf ). The ratio is accurately unity with- 
out dependence on declination, implying that the PSF modeling is 
accurate. 
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Fig. 6. — Star/galaxy separation based on PSF— model magni- 
tude, r pa f — r mo[ j i. A model magnitude is the best galaxy deVau- 
couleurs or exponential profile convolved with the PSF in the sense 
of describing the data. In the event of a mismatch of the modeled 
PSF with the real PSF, or extended light, the model magnitude 
will measure more light than the PSF. It is thus a sensitive indi- 
cator of whether an object is well fit by the modeled PSF. Objects 
classified as stars are marked in red, galaxies in blue. There is a 



clear stellar locus at r 



psf 



odel ~ 0- That the stellar locus bends 



towards negative values at r pa f > 22 is likely due to errors in the 
PSF fit. 



1 + exp((m — jU)/<r) 

where [i is the magnitude limit of the catalog (defined to 
be the magnitude at which C{m) = 50%), /o is a normal- 
ization constant and the parameter a controls how fast 
the completeness falls when it reaches the completeness 
threshold. We use this fitting function to determine the 
depth of the coadd galaxy and star catalogs, /j,g and \is 
respectively. 

Our results are illustrated in Fig. [7] The plots on the 
first row show the r-band completeness for the same two 
fields pictured in Fig. 2} On the right we have the coadd 
field (run 206, camcol 3, field 505) and for comparison, 
on the left, one of the 28 single pass images used as input 
for that particular field. The coadd reaches r — 24.3 for 
point sources and r = 23.4 for galaxies, going about 2 
magnitudes deeper than a single pass image, as expected. 
The plots on the second and third rows show the coadd 
results for the other filters. Table [5] is a compilation 
of the median and rms values of the coadd magnitude 
limits, calculated for 50 fields randomly selected. 

5.5. Color-Color Diagrams 

Stars populate a well-defined locus in the color-color 
space almost independent of magnitude. Therefore color- 
color diagrams are useful to access the quality of the 
photometry. Figure [8] shows the color-color diagrams of 
a isolated and well measured sample of coadd stars se- 
lected at —5° < RA < 0°, in a high galactic latitude 
in Stripe 82 (see query in the Appendix). The sample 
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Fig. 7. — Completeness as function of r-band model magnitude for the coadd run 206, camcol 3, field 505 (corresponding to the image 
shown in Fig. [21. Point sources are represented in red, galaxies in blue. The solid lines are our measurements and the rms uncertainties 
are representea as light-colored regions. Dashed lines are the best fit model. There is a ft: 0.09 mag scatter in the measurement of the 
completeness level as measured from frame to frame. The numbers refer to the 50% completeness level which we find to be a much more 
stable fit than the 95% level. Top: Comparison between the coadd (right) and one of the single pass images used in the coaddition (left), 
showing that we achieve ~ 2 mag deeper as N = 28 images contributed into this particular frame. Middle and bottom: Results for griz on 
the same frame, showing the typical depth of the coadd in each bandpass. 
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TABLE 5 

50% Completeness Limits in the 
Coadd 



filter 


MS 




MG 


"■(mg) 


u 


23.63 


0.06 


23.25 


0.23 


g 


24.56 


0.10 


23.51 


0.18 


r 


24.23 


0.08 


23.26 


0.14 


i 


23.74 


0.15 


22.69 


0.17 


z 


22.29 


0.09 


21.27 


0.23 



Note. — Coadd magnitude limits for 
stars (fJ-s) an d galaxies (hg)- Values re- 
ported are medians and rms calculated 
for 50 fields randomly selected accross 
Stripe 82. 



was split into 1 magnitude bins in r-band. At brighter 
magnitudes the intrisic thinness of the stellar locus is 
apparent. At fainter magnitudes statistical noise begins 
to dominate. Since this is an r-band sample, the it-band 
stars in the 20 < r < 21 panel are quite faint; one cannot 
read from these diagrams where the signal to noise of the 
data degrades. One can read in these diagrams just how 
good the photometry can be for clean samples of stars. 

For a reasons which we have not unearthed, the sat- 
urated flag did not propogate through PHOTO. Beware of 
objects with m < 15.5. 

5.6. Number Counts 

The number counts of stars and galaxies allow us to 
asses the depth of the coadd and the success of the 
star/galaxy separation. Figure [9] shows the «-band star 
and galaxy counts in regions of 5 deg 2 patches at a variety 
of galactic longitudes and latitudes along Stripe 82. The 
galaxy number counts show the Euclidean m 6 power 
law expected at i < 20, the slow change of slope due 
to cosmological volume and galaxy evolution at fainter 
magnitudes, and a roll off at i « 23.5 due to completeness 
issues. We have seen in panel 5 of Fig. [7] that the galaxy 
catalog is ~ 50% complete at i — 23, consistent with the 
deviation from the slowly rolling power law index seen 
here. We conclude that this plot shows nothing seriously 
wrong with the galaxy counts from 16 < i < 23.5. 

The star counts are going to be slightly more problem- 
atic. They are rough power laws that cross the galaxy 
counts at i ~ 20. The expected roll-over from incom- 
pleteness is apparent near i w 23.5. The steep rise in 
counts of stars seen easiest in the RA = +57° patch data 
points is from galaxies being misclassified as stars: the 
galaxies outnumber the stars by > 30 so even a small 
misclassification rate result in a large number of galaxies 
scattered into the star class, and one can see in Fig.[6]that 
at faint magnitudes the galaxy locus crosses the stellar 
locus. Once it does, stars can no longer be distinguished 
by galaxies using psf-model alone. 



The models are from Trilegal star count modellingj Gi- 
rardi et al. 2005 online version vl.4). In F igure 9| the 
galact ic model parameters are those found by |Juric et al. 
( 2008 ) using their SDSS star count tomographic mapping 
technique. The slope of the rough power-law is due to 
a combination of thin disk at brighter magnitudes and 
halo stars at fainter magnitudes. This model doesn't 
fit the data particularly well at the RA = —31° patch, 
or at lower RA. Models with a lower exponential scale 



height, such as those derived from SDSS data using Tri- 
legal model fits by B. Santiago ( private communication' 
or from the SDSS m-dwarf fits of Bochanski et al. ( 2007 ' 
do describe the Stripe 82 data at RA < -31°. Figure[Tc 
shows number counts in all five filters for the RA = — 31 5 
field along with bands showing the range of star counts 
from the Santiago and Bochanski galactic models. This 
figure also shows, in all filters, the spike in counts as 
galaxies enter the star catalog in large numbers at faint 
magnitudes, and often a flattening of star counts at mag 
~ 22. As we have seen in section §5. 3| stars are entering 
the galaxy catalog in higher than expected numbers at 
these magnitudes due, most likely, to a PSF modelling 
problem, and the flattening is probably a reflection of 
this. Nevertheless, the data seems to indicate that there 
is an interesting problem in combining the deep Stripe 82 
star counts with the single pass SDSS star counts. Over- 
all, this and the previous figure show reasonable agree- 
ment between the Stripe 82 star counts and the mod- 
els. Comparison with star count models in this paper 
is an approximate attempt to validate the overall sanity 
of t he stellar counts. F or a a detailed Galaxy study, see 
e 
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6. SCIENCE 

The coadd galaxy catalog is the largest homogeneous 
photometric catalog of its depth. It has 13 million galax- 
ies and we have shown that it is complete to r = 23.5, 
being 2 magnitudes deeper than the SDSS single scan 
imaging survey. This implies that, despite covering an 
area 40 times smaller, the coadd volume has roughly 1/3 
of the total volume of the SDSS data on the Northern 
Galactic Cap. The scientific questions that can be ad- 
dressed by exploring this wealth of data are numerous 
and diverse. In this section we present a brief overview 
of recent and ongoing analyses using the coadd data. 

The presence of deep imaging provides a template for 
variable object measurements and a means to check on 
the fainter objects of the si ngle scan SPSS data. The 
SDSS-II Supernova Survey (Frieman et al. [2008) used 
the coadd images as a high-quality, photometrically cal- 
ibrated template for carrying out image subtraction to 
discover supernovae. Later observations of the host 



galaxy often start with t he coadd images ( Lampeitl et al. 
2010 Gupta et al.||2011 ) to locate the host. As much of 
the stripe 82 data was taken as part of the SDSS-II SN 
survey, this use is expected. 

There are many uses that are less tightl y related to the 
goals of taking the data. For example, Liu et al. ( 2011[ ) 
use the greater depth of the coadd as a means to check 
their science at the main survey limits, as they search for 
paired quasars in the SDSS and use the coadd to check 
for the fraction of interacting pairs missed due to the 
surfac e brightness limit of the single scan. Jiang et al.| 
(2008) used the coadd to perform the the discovery of 
5 quasars at redshift ~ 6 and t o compute th e high red- 
shift quasar luminosity function. Vidrih et al. ( 2007 ) find 
white dwarf candidates. Our exploration of the Trilegal 
models in section §5.6| shows that these data could be 
used to constrain Galactic models. 

Our main interest is in cosmology, and here there is a 
clear path of work. Photometric redshift measure ments 
for the coad d galaxies were obtained by our group (Reis 



et al. 2011) using a neural network trained on spectro- 
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Fig. 8. — Color-color diagrams of stars in a 12 degree 2 field, in bins of 1 magnitude from 18 < r < 24. Photon noise begins to dominate 
the bluer color in the u—g vs g — r plot at r > 20, the g — r vs r — i plot at r > 21, and the z band in the i — z vs r — i plot at r > 21. 
As the stellar locus is intriniscly very thin, these plots provide a sense of how good the photometry is at brighter magnitudes, and where 
statistical noise begins to dominate. 
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Fig. 9. — Star and galaxy counts in 5deg 2 patches: the data are 
points with poisson errors, and models are shown by lines. The 
three patches are at RA = -31°, Dec = 0° (I = 58°, b = -40°), 
RA = 13°, Dec = 0° (I = 123°, b = -63°), and RA = 57°, Dec 
= 0° (I = 188°, b = -40°). The galaxy counts (in black) show the 
expected Euclidean power law at i < 20, the slow change of slope 
due to cosmological volume and galaxy evolution at 20 < i < 23.5, 
and the roll off at i m 23.5 due to completeness issues. The star 
count models fit the data well enough for the desired purpose, 
until i w 23.5, where there is a sudden upturn in star counts, most 
evident in the RA = 57° data. This is certainly due to galaxies 
being classified as stars. 

scopic redshifts obtained on Stripe 82. As Stripe 82 
is easily accesible by telescopes in both the North and 
South hemispheres, t he area has been well studied spec- 
troscopically. For the |Reis et al. ( 2011 ) work, we used the 
Canadian Network for ( J bservational C osmology Field 
Galaxy Survey (CNOC; |Yee et al.||2000| ), the Deep Ex- 



traga lactic Evolutionary Probe (DEEP2; |Weiner et al. l 
2005|), the WiggleZ Dark Energy Survey (Drinkwatcr 



et al. 12010 ), and the Visible imaging Multi-dbject Spec 
trographvery La rge Telescope Deep Survey (VVDS; Le 



Fevre et al.|2005|. The mean photo-z error acheived was 
a(z) = O.OiJl and the photo-z catalog is reliable out to 
- 0.8, after which the low S/N of obj ects in the SPSS 
■band becomes the limiting factor (see; Reis et al.|2011 



for a detailed discussion) . This photo-z catalog has been 
made public at the same time as this paper as a value- 
added catalog. 

Since the photometric redshifts go to z 0.75, an in- 
teresting program is cluster finding in the range 0.5 < 
z < 0.75. Preliminary cluster catalogs have been pursued 
by our group using the Gaussian Mixtu re Brightest Clus- 
ter Galaxy (GMBCG; |Hao et al.|[20To| ) and the Voronoi 



Tessellation clust er finder in 2+1 dimensions ( Soares 
Santos et al.|2011 VTT) algorithms. We plan on pursing 



a search for blue clusters using the VTT as a finder and 
the GMMBCG as a red sequence measuring engine aim- 
ing at studying clust er formation and evolution. 

In Lin et al. (2011 ) we report the measurement of cos- 
mological parameters from the cosmic shear signal in the 
coadd. All weak lensing analyses require accurate shape 
estimation parameters, but the cosmic shear is an ex- 
treme case, due to its very low signal. PHOTO measures 
second moments and related parameters needed for weak 
lensing, but several systematic errors in the PSF had to 
be corrected. Some of these systematic were doubt due to 
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Fig. 10. — The Trilegal model star counts in u,g,r,i,z and coadd 
star counts at RA = -31°, Dec = 0° (b = -40°, I = 58°). The 
bands show a reasonable range of models, from one derived from 
Bochanski, the other from a private communication model of San- 
tiago, modified to have the thin disk scale height z = 80 pc. 

our mis-modelling of the PSF because of the mis-matched 
weights (see section \ 5.2 ). After the corrections and weak 
lcnsing-spccific quality cuts, the coadd data provides ~ 6 
galaxies per arcmin 2 for the analysis, 6 times more than 



the SDSS single pass data. Huff et al. (2011), in an in- 
dependent work, have also measured the cosmic shear on 
Stripe 82. They did not use the coadd described in this 
paper, but instead, made their own coadd optimizing for 
weak lensing. 

A related program is to measure the masses of clusters 
found in th e SDSS Coadd area. The MaxBCG ( |Koester| 
et al.|[20"07|) cluster catalog overlaps with Stripe 8 2 using 



single pass d ata (as do oth e r cluster catalogs, e .g., Geach 



et al.| ( |2011[ ); |Dong| ( |2011| ); [S"zabo et al.| ( |2011[ >). 

We have performed a stacked cluster weak lensing anal- 
ysis with the MaxB CG clusters as lens es and the coadd 
galaxies as sources (Simet et al. 20111. We divide our 
cluster sample in bins of richness and measure a mass- 



richness rel ation consistent with previous work ( Johnston 
et al.|2007 ). This demonstrates that the coaddititon pro- 
cess does not dilute the lensing signal. As we detect an 
increasing sig nal as a fun c tion o f source redshift we also 
conclude in |Simet et al. ( |2011| that 

signal ii 



weak lensing tomography signal in the coadd 



have detected 
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Since the survey covers a large area on a part of the 
sky that has been heavily studied there are also many 
opportuniti es for multi-waveleng th studies, from the x- 
ray (XMM: |Mehrtens et al.|2011 ) to the m i crowave (SZ: 



Menanteau et a 



2011||Sehgal et al. 



20101 JHand et al.|[20TT| |Reese et al.] 



2011) 



Ttis unlikely that we have surveyed all of the science 
the coadd has been put to use for, as it is part of the 
DR7 database and easily accesible by the community. 
This paper serves as as a technical description of this 
widely available dataset. 

7. CONCLUSIONS 

The SDSS performed repeat scanning of the equato- 
rial region in the South Galactic Cap known as Stripe 
82. The amount of data observed was comparable to the 
single scan coverage of the SDSS footprint. The aim of 
the work described here was to coadd this data and an- 
alyze it using the SDSS pipeline framework, notably the 
PHOTO pipeline. Roughly a third of the existing Stripe 
82 data was coadded, limited in time by when the work 
was performed. The runs included calibrated and uncal- 
ibrated data, so a relative calibration scheme was devel- 
oped and applied. The images were mapped onto a SDSS 
run format output grid using the SDSS astrometry, were 
coadded using a S/N weighting that includes seeing, and 
inverse variance maps computed. The coadded images 
were run through PHOTO using PSF models computed by 
coadding the PSFs. 

The resulting catalogs have median seeing in r of 1.1" 
and varies band to band following Kolmogorov scalings. 
The catalogs are 50% complete to r = 23.5 (galaxies) 
and r = 24.3 (stars). The photometry is good to .5% in 
g,r,i, and 1% in u and z, as measured against the llvezic] 
et al. (2007) star catalog. The PSF is modeled and de- 



spite minor issues it is useful for precision photometry. 
Color-color diagrams of stars show a sharp and thin stel- 
lar locus, and the number counts of stars at a variety of 
positions show agreement with reasonable galactic mod- 
els. There are identifiable regions in psf-model vs magni- 
tude space where stars are being misclassified as galaxies, 
and we give suggestions on how to eliminate these. Thus 



we have constructed high quality SDSS catalogs on the 
Stripe 82 region from images that are two magnitudes 
deeper than the single pass SDSS data. 

The coadd catalog is largest homogeneous precision 
photometric catalog complete to r=23.5. The catalog 
has 13 million galaxies, and has a variety of uses from 
galactic structure to large scale structure and weak lens- 
ing, to cosmology. The data, including both images and 
catalogs, are available through the standard SDSS dis- 
tribution channels. 
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APPENDIX 

QUERY FOR CLEAN PHOTOMETRY STARS 

Here we provide the query used to obtain the sample of isolated and well measured stars used for the color-color 
diagrams. The query is to be run on the Stripe 82 database of the SDSS Catalog Arquive Server (CAS). 



SELECT 

ra, dec, run, camcol, field, 
u, g, r, i, z, 

psfMag_u, psfMag_g, psfMag_r, psfMag_i, psfMag_z, flags, 

psf magerr_u, psfmagerr_g, psfmagerr_r, psf magerr_i , psfmagerr_z 

FROM 
PhotoObjAll 

WHERE 

((flags & 0x10000000) != 0) 
AND ((flags & 0x8100000c00a4) = 0) 
AND (((flags & 0x400000000000) = 0) or 

(psf magerr_r <= 0.2 and psf magerr_i<= 0.2 and psfmagerr_g<=0 . 2) ) 
AND (((flags & 0x100000000000) = 0) or (flags & 0x1000) = 0) 
AND (run = 106 or run = 206) 
AND type = 6 
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AND mode = 1 

AND ra between 355 and 
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